Signal analysis by expansion over the squared eigenfunctions of an associated Schrodinger 

operator 
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Abstract 

This article introduces a new signal analysis 
method. The main idea consists in interpreting a 
pulse-shaped signal, after multiplying it by a posi- 
tive parameter, potential of a Schrodinger op- 
erator and representing this signal with the discrete 
spectrum of this operator. We present some results 
obtained in the analysis of the arterial blood pressure 
with this method. 

Introduction 

Let H(V) be a Schrodinger operator in L 2 (R): 



d 2 

H(V)%b = -^V' + Vif>, ip G V(H(V)) = H 2 

where V, called a potential, satisfies: 

am 

V G L\(R), 
with 

L\(R) = {V\ 



dx r 



(1) 



UgL^M), to =1,2, (2) 



|V(s)|(l + \x\)dx < oo}. (3) 



H 2 (R) refers to the two order Sobolev space. The 
spectral problem of H(V) is given by: 

d 2 V> 



dx 2 



+ V(x,t)i/> = k z ip, keC , ieR, (4) 



where k 2 and ip are respectively the eigenvalues of 
H(V) and the associate eigenfunctions. Under hy- 
pothesis (j2j), H(V) is self adjoint and its spectrum 
consists of: 

• a continuous spectrum equal to [0, +00), 

• a discrete spectrum composed of negative eigen- 
values of multiplicity 1. 

Let y be a real valued function representing the 
signal to be analyzed such that: 



y G L\(R), y(x) > 0, Vi£ 
d m y 



then H(—xy) is defined in L 2 (IR) for all x > by: 

H(-xv) = -^2-XV, V(H(- X y)) = H 2 (R). (6) 
dx L 

Under hypothesis (|5|), there is a non-zero (PQ, 
corollary 2.4.4), finite number N x ([2], theorem 1) 
of negative eigenvalues of the operator H(~xy)- 

1 Semi-classical signal analysis (SCSA) 

We define the SCSA approximation by: 

Definition 1.1. Let y be a real valued function satis- 
fying hypothesis f5|) and x a, positive parameter, then 
the SCSA approximation is defined by: 



— K «xV'nx( X )' X ^ 

X 



(7) 



n=l 



where i/ie negative eigenvalues of H{~xy) 

with n nx > and ki x > K2 
1, • • • , iV x and ipnx> n = 1, - ■ ■ 
L 2 -normalized eigenfunctions. 



, A^ x are i/te associate 



The values — — can be interpreted as particular 

X 

values of a signal that can not be extracted using 
usual sampling methods: values on the discretisa- 
tion points and extremal values. The semi-classical 
interpretation related to the Bohr-Sommerfeld quan- 
tification seems to be natural in this case [3J. Let us 
illustrate this idea briefly. Suppose that y, Vx G M. 
is a positive signal that vanishes rapidly. We denote 

1 k 2 
h = and = — . Then the spectral problem of 

Vx " x 

H{~xy) can be written in the following form: 



dx r 



1,2, 



(5) 



-^ 2 -^(\AftX x)-y(x)7p(x/\i~x, x) = X h ifj(x/\i~x, x). 

(8) 

When x ~~ * +00, hence h — > 0, the problem is equiv- 
alent to a semi-classical problem [4j. If < y(x) < 
Umax, Vx G R, then it is well-known that the negative 




eigenvalues of the operator —h 2 ——^ — y(x), namely 

dx z 

K 2 
72V 

, n = 1, • • • , N y are comprised between —y ma x 

X 

and as it is illustrated in figure [TJ For a fixed value 
of x, they correspond to particular values of —y. 
Therefore, we can associate to a signal, taken 

potentiel well, the values — — , n = 1, • • • ,N y which 

X 

are used as some representative values of the sig- 
nal. This is a new quantification approach: a semi- 
classical quantification. It is based on some standard 
results from semi-classical analysis gathered in the 
following proposition (for more details, see [3]): 

Proposition 1.1. i) Let y be a function sat- 
isfying hypothesis |5|J. Let % > and —K^x' 
n = 1, ••• ,N X with — «i x < ~k% x < ••• < 
the negative eigenvalues of XV)- We sup- 
pose that < y{x) < y ma x, Vx G M then, 

K 2 

— Umax i — 1) ' ' ' ; N y . 

x 

ii) Moreover if y £ C°°(R) such that for one 
7o € R, miniR (— y + 70) > and for all a G N, 
there is a constant C a > swc/i £/iaf < 
C a (— y + 7o) ; ^en ei;ery regular value of y is an 

K 2 

accumulation point of the set f-^, X > 0; n = 
1, • • • , N x ) (v is a regular value if < v < y ma x 
and ify(x) = v then \ dy j^ \ > 0). 

2 Numerical results 

Figure [2] shows the reconstruction of one beat of 
an Arterial Blood Pressure (ABP) signal with the 
SCSA. We noticed that only 5 to 10 negative eigen- 
values are sufficient for a good reconstruction of an 
ABP signal. The SCSA was also applied for the 



separation of the systolic and the diastolic pressures 
which describe fast and slow phenomena respectively 
[5]. We also point out that the SCSA introduces 
some interesting parameters that give relevant phys- 
iological information. These parameters are the neg- 
ative eigenvalues and the so called invariants that 
consist in some momentums of n nx , n = 1, • • • , N x 
[6]. For example, these new cardiovascular indices al- 
low the discrimination between healthy patients and 
heart failure subjects [6|. 
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